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RECOVERY OF STOKES WAVES FROM VELOCITY MEASUREMENTS ON 

AN AXIS OF SYMMETRY 

BOGDAN-VASILE MATIOC 


Abstract. We provide a new method to recover the profile of Stokes waves, and more generally of 
waves with smooth vorticity, from measurements of the horizontal velocity component on a vertical 
axis of symmetry of the wave surface. Although we consider periodic waves only, the extension to 
solitary waves is straightforward. 


1. Introduction 

A Stokes wave is a two-dimensional periodic wave that travels at constant speed over a flat 
bed solely under the influence of gravity. The flow beneath a Stokes wave is irrotational and the 
surface profile is symmetric with respect to both crest and trough lines. Initiated by Stokes [30], 
the rigorous study of Stokes waves was developed steadily to uncover many fascinating aspects of 
these remarkable water flows, cf. e.g. [2, 6, 8, 11, 33]. 

The problem of reconstructing the free surface of regular Stokes waves from measured flow data, 
such as the pressure on the fluid bed, is of great importance in the field of fluid dynamics, cf. 
[1, 5, 13, 14, 24, 28, 34]. A very simple relation between the pressure on the bed and the elevation 
of the water surface is obtained when assuming that the pressure is hydrostatic, see the discussion 
in [5, 28]. This method is used for tsunami detection, but, as it does not account for nonlinear wave 
effects, prediction errors are quite frequent [14]. On the other hand, the transfer function approach 
[16, 23] uses a linear formula relating the wave surface elevation to the pressure on the fluid bed, 
a large body of research being dedicated to determining a suitable transfer function. This latter 
method restricts primarily to linear irrotational waves. Therefore, it does not account for nonlinear 
effects and sometimes overestimates the height of large waves [34]. 

Within the exact theory of irrotational water waves, a nonlinear and nonlocal relation between 
the wave profile of solitary waves and the pressure on the bed was derived in [13, 28]. An alternative 
exact local recovery formula was obtained in [7, 20]. For a recovery formula in the setting of steady 
periodic waves we refer to [4, 5]. In fact it is theoretically possible to recover the profile of solitary 
waves from the pressure on the bed even when allowing for rotational flows, provided that the 
vorticity distribution is real-analytic [19]. 

The goal of the present note is to present a new exact method to recover the profile of periodic 
traveling waves and the velocity field within the flow from measurements of the horizontal velocity 
on a vertical axis of symmetry of the wave. Our analysis is within the setting of rotational periodic 
waves with a smooth vorticity distribution, the extension to solitary waves is straightforward. In 
particular, our method can be used to reconstruct the free surface of Stokes waves. That the profile 
of rotational periodic gravity waves with a smooth vorticity is determined by the horizontal velocity 
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on a vertical axis of symmetry - which is either a crest or a trough line - is quite intriguing and 
it underlines the strong connection between the flow and the free surface motion. This intricate 
connection is suggested also by the result of [26] where it is shown that the symmetry of the free 
wave surface can be characterized intrinsically in terms of the underlying flow. 

The main tool used in our approach is a very recent observation pertaining to the regularity of 
rotational flows: any weak solution of the water wave problem possesses real-analytic streamlines, 
the radius of analyticity being independent of the position in the fluid, cf. [15]. Finally, we em¬ 
phasize that the horizontal velocity on crest or trough lines has been measured in many laboratory 
experiments [18, 21, 22, 29, 31] and therefore we expect that the practical accuracy of our recovery 
method can be verified numerically. 


2. Equivalent formulations of the water wave problem 


We consider steady periodic water waves propagating at the surface of a two-dimensional fluid 
body with wave speed c > 0. The Cartesian coordinates (x, y, z) are fixed such that the x-axis is 
the direction of wave propagation, the y -axis points vertically upwards, and the flow is independent 
of the ^-coordinate. The fluid layer is bounded from below by the flat bed y = —d, whereby d > 0. 
Neglecting viscosity effects and temperature variations, the motion of the fluid is governed by the 
Euler equations. As we consider steady waves only, the space-time dependence of the free surface, 
of the pressure, and of the velocity field has the form [x — ct). Eliminating time by the change of 
coordinates (x — ct,y,z ) H > ( x,y,z ), in the new coordinates system in which the origin moves in 
the direction of propagation of the wave with wave speed c, the wave is stationary and the flow is 
steady, the equations of motion being encompassed in the following system 


( 2 . 1 ) 


Hereby, y = rj(x) is the surface of the wave, := {(x,y) : —d < y < 77 ( 2 )} is the fluid domain, P 
denotes the dynamic pressure, u is the horizontal velocity, v is the vertical velocity, g is the gravity 
of Earth, the constant Po is the uniform air pressure, and the water’s density is set to be p = 1. 
Classical solutions ( u,v,P,r /) of (2.1) belong to the following regularity class 


c)u x + VUy 

= -p x 

in 

c)v x + VVy 

II 

1 

1 

in fijj, 

U x + Vy 

= 0 

in 

P 

= Po 

on y = rj(x) 

V 
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1 

II 

on y = rj(x) 

V 

= 0 

on y = —d. 


(ix, v, P, rj) € (C^ er (Q v )) 3 x C } 


'per 


( 2 . 2 ) 


the subscript per being used to express the A—periodicity of (u, v, P, rf) with respect to the variable 
x , and where the constant A > 0 is the wavelength. A last assumption we make is that the wave 
speed c exceeds the horizontal velocity of each particle within the fluid 


sup u < c. 

Q77 


(2.3) 


The formulation (2.1) can be also recast in terms of the streamfunction ^ which is defined by 

/ y _ 

(■ u(x,s)-c)ds for (x,y) € £l v , (2.4) 

-d 
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the constant po < 0 , which represents the relative mass flux, being given by 



(2.5) 


It follows from the third and the last two equations in (2.1) that po is independent of x. Let us 
also observe that since \7ip = (—v,u — c), the vorticity of the flow is given by the scalar function 
co := u y — v x = Aip. 

Because of (2.3), the map T~L : 12^ —> 12, with 12 := M x (p o, 0), defined by the relation 


H(x,y) ■= {q,p)(x,y) = (x,-ip(x,y)) for (x,y)en v ( 2 . 6 ) 

is a diffeomorphism, that is Ti G Diff 2 (12 J? ,12). Moreover, it follows from (2.1), cf. e.g. [25, Section 
2], that cj o =: 7 depends only on the variable p. The function 7 G C([po>0]) is the s °~ 
called vorticity function. Note that both cj and 7 are zero constant functions for Stokes flows. 
Using Bernoulli’s law, we find that (p,ip) £ Cp er (®0 x Cperi^v) solves the following free boundary 
problem 

Aip = 7 {—ip) 

IW >| 2 + 2g(y + d) = Q 
ip = 0 
ip = -po 

with the positive constant Q := 2 (E — Pq). In terms of ip the condition (2.3) appears as 


in fijj, 

on y = r)(x), 
on y = r](x), 
on y = —d, 


(2.7) 


swpip y <0. ( 2 . 8 ) 

^77 


For the existence theory, but also when studying the properties of solutions of (2.7), cf. [12, 15], 
it is useful to express the hydrodynamical problem in terms of the height function h : Q —y R that is 
defined by h := yoPT 1 +d. The function h G C' 2 er (n) associates to each pair (q,p) G O, the height of 
the water particle (x,y) := 'H~ 1 (q,p) above the bed. Particularly, we have that rj(x) = h(x, 0) — d 
for all x G M.. It follows easily from the definition of h and (2.7), cf. e.g. [12], that h solves a 
quasilinear elliptic equation with nonlinear boundary conditions 

(1 T hq)hpp 2 hphqhpq + hphqq 7 hp — 0 in n, 

( 1 + hq + (2 gh — Q)hp = 0 on p = 0, (2-9) 

h = 0 onp = pq, 

while the assumption (2.3) is now equivalent to 


inf hr, > 0 . 

n 


( 2 . 10 ) 


In the next Section, we use the equivalence of the three formulations within the set of classical 
solutions. The proof of this result follows easily from the arguments presented in [ 6 , Chapter 3] and 
in [25]. 


Proposition 2.1 (Equivalence of the three formulations). The following are equivalent: 
(i) the equations (2.1) and (2.3) with ( u,v,P,rj ) satisfying (2.2); 

(ii) the equations (2.7) and (2.8) with (r),ip) G Cp er (M) x Cp er (Q v ) and 7 G C([po,0]); 
(Hi) the equations (2.9) and (2.10) with h G Cp er (Q) and 7 G C([po>0]). 
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3. The main result and discussion 

In order to state our main result, we emphasize first that any (vertical) axis of symmetry of the 
wave must either be a trough or a crest line. This is a simple consequence of the real-analyticity of 
the wave profile, a property which holds a priori for solutions in the class (2.2), cf. [15, Corollary 
1.2]. The main result of this paper is the following theorem. 

Theorem 3.1. Consider a regular 1 two-dimensional symmetric and periodic gravity wave with 
vorticity ui G C 00 (12^) that travels over the flat bed y = —d with wave speed c > 0. Furthermore, 
assume that the horizontal component u of the velocity field is known on the line of symmetry 

{(x 0 ,y) : —d < y < y{x 0 )} 

of the wave. Then, we can recover the wave profile, the velocity field and the pressure distribution 
within the flow. 

Before proving Theorem 3.1 it is worthwhile to add some remarks. 

Remark 3.2. (a) Theorem 3.1 is true in particular for Stokes waves for which the vorticity is 

zero. We emphasize that a constant nonzero vorticity, for which the assumption u G C'°°(fl T? ) 
is also trivially satisfied, is relevant in many physical situations, cf. the discussion in [32]. 

( 6 ) We remark that to the best of our knowledge, there is no rigorous proof for the existence of 
two-dimensional periodic steady gravity waves which are not symmetric. For criteria which 
ensure the symmetry of such waves with respect to both trough and crest lines we refer to 
[9, 10, 26, 27], 

(c) In the proof of Theorem 3.1, we do not use of the assumption of periodicity. Particularly, 
our result shows that also the surface of solitary waves can be recovered when knowing the 
velocity field on a axis of symmetry. 

( d ) The values of u on the vertical symmetry line [x = xo] and the wave speed are not both 
needed, but just the values of u — c on this axis of symmetry. 

(e) The proof of Theorem 3.1 uses to a large extent the fact that the streamlines of the flow 
are all real-analytic, the radius of analyticity being uniform in the fluid. An open question 
is whether the radius of analyticity can be chosen to be oo. 

Proof of Theorem 3.1. Let 77 , ( u , v), and P be the unknown wave profile, velocity field, and pressure 
distribution, respectively. Because of the equivalence result in Proposition 2.1, we only need to 
determine the height function h. Observing that the problem (2.1) is invariant with respect to 
horizontal translations we may set xq = 0. Hence, y is symmetric with respect to the vertical line 
[x = 0], 

Our assumptions together with (2.5) enable us to determine the relative mass flux constant po, 
cf. (2.5). Knowing the strip where the height function h G Cp er (fl) is defined, we next show that 
h is well-determined by the restriction of (u — c ) on the axis of symmetry [x = 0]. To this end, 
we recall that Aip = uj in and that if is constant on dd v . From our assumption uj G C'°°(H^) 
and the real-analyticity of the wave profile we deduce that if G C'°°(n r? ), cf. [17, Theorem 6.19]. 
Therewith, we get 7 = uj o 'H~ x g C°°([po, 0]) and also h G C£f r (Q). 

^By regular wave we mean a solution (u,v,P,g) of the hydrodynamical problem (2.1) that satisfies (2.2) and (2.3). 
Notice that this class of solutions ensures that the three formulations are equivalent, cf. Proposition 2.1. The 
condition (2.3) is satisfied by regular Stokes waves automatically, only the Stokes wave of greatest height having the 
property that u = c at the crest where the free water surface forms an angle of 120 degrees, cf. [33]. 



RECOVERY OF STOKES WAVES FROM VELOCITY MEASUREMENTS 


5 


Using the weak elliptic maximum principle in the context of (2.9) as in [26, Lemma 3.2] shows 
that the symmetry of the free water surface with respect to [x = 0], that is r]{x) = p(—x) for all 
x £ R, also implies the symmetry of h, that is h(q,p ) = h(—q,p) for all (q,p) £ U. Even more holds: 
there exists a constant L > 0 such that 


\\d™h\\ c2m < L m l (m — 2)! for all m> 2, (3.1) 

cf. [3, Theorem 3.3] (see also [15] for a more general result). The estimate (3.1) ensures the real- 
analyticity of the streamlines as [q i—>• h(q,p)\, p £ [po,0], is real-analytic. Indeed, (3.1) implies for 
each (qo,p) £ that 


h(q,p) ~ Y 

k=0 


dgh(q 0 ,p) 

k\ 


(q-qo) k 


/ Iloilo „ t n+i^(L\q-qo\r +1 
-^TW lq ~ qo1 - - L - 


(3.2) 


if L\q — qo\ < 1. In particular, the radius of convergence of the Taylor series is independent of the 
point ( qo,p ). This latter property enables us to extend the function h to the whole strip U, provided 
that we know it on the subset [— (2L) _1 < q < [2L)~ l ] C Q, where the Taylor series of [q i-A h(q,p)\ 
converges for each p £ [po,0]. 

With this, our task reduces to determining h in the rectangle [—(2L)” 1 < q < (2L) -1 ]. We note 
that the symmetry of h implies in particular that dg n+1 h(0,p) = 0 for all n € N and p £ [po,0]. 
This property and (3.2) yield that 


h(q,p) = Y a2n ^ q2n for M < ( 2L ) 1 and V G [po,0], (3.3) 

n=0 

whereby a^n '■= ((2n)!) _1 0g n /i(O, •) £ C°° ([po, 0]) for all n £ N. Similarly to (3.2), one can deduce 
from (3.1) that 


OO OO CO 

h p {q,P ) = Y a 2 n(p)q 2n i h q (q,p) = Y 2 na ‘ 2 n(p)q 2n ~ 1 , hp P (q,p) = Y a * nW? 2 ", 

n =0 n =1 72=0 

CO CO 

hqp(q,p ) = Y 2 na 2n(p)q 2n ~ 1 ^ h qq (q,p ) = ^ 2n(2n - 1 )a 2 n (p)q 2n ~ 2 

n =1 n =1 


(3.4) 


for |q — qo\ < (2 L)~ l and all p £ [po>0[. 

Since the values of (u — c ) on the symmetry line \x = 0] are known, we can determine the 
coefficient function ao as follows. From the definition of h we obtain that h{ 0, —ip(0,y)) = y + d for 
all y £ [— d, 77 ( 0 )], whereby, in virtue of (2.4) and (2.5), 

V>(0, y) = ~Po + / (n(0,s) -c)ds, y£[-d,r]( 0)]. 

J-d 


Because of (2.8), the map — r/;(0, •) : [~d,rj( 0)] [po 5 0] is bijective and its inverse is thus also 

known. With the notation from Section 2 it is the function y o %~ 1 (0, •) : [po,0] —> [—d,r]( 0)]. 
Hence, ao is determined by the formula ao = y o ?^” 1 (0, •) + d : \po, 0] —>• [0, d + 77 ( 0 )]. 

To finish the proof, we show that ao determines all the other coefficient functions 02 ^, n £ N, 
n > 1. Indeed, plugging the Taylor expansions (3.2) and (3.3) into the first equation of (2.9) yields, 
after identifying for a fixed p £ [po; 0 ] the coefficients of each power of q , that 


Aq + Co + Dq — 0 and A 2 n T B 2 n T C 2 n + D 2 n — 0 for n > 1 . 


(3.5) 
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Hereby, we used the following relations 

oo oo oo 

(1 + h 2 q )h pp = ^2 Mnq 2n , - 2 h q h p h pq = B 2n , Q 2n h p h qq = ^2 C 2n q 2n , ~ 7 ^p = 22 D2nQ 

n —0 n =1 72=0 

with 

n k 

Mn = a'ifi + 22 a 2 (n-k) ( 22 - l + l)a 2 za 2 (fc-Z+i)) , 


2 n 


n =0 


k=l 
n—1 


1=1 
k+1 


B 2 n — 2 ^ a 2fa-fc-i'i ( 4£(fc / + 2)a 2 za 2 (fc-z+2) 

fc =0 Z =1 

n n—k 

c 2n = + l)( 2 fc + 2 )a 2 fc +2 f a 2l a 2(n-k-l) 


k =0 


Z=0 


n—k 


D 2 n — -1^2 a ' 2k ( a U a 2(n-k-l)) ■ 
k =0 Z=0 

The first equation of (3.5) is equivalent to 


«2 = 


7°o «o 


2 ( oq ) 


/ t 2 ’ 


(3.6) 


whereby, in view of u < c in , we have Oq > 0. Moreover, from the second equation of (3.5) we 
find that 


1 


fl 2 n +2 — 


(2n + l)(2n + 2)(a' 0 ) 2 


~ a 2n ~ 


22 fl 2 (n-k) ( 22 4 l{k - l + l)a 2 za 2 (fc_; + i) 


k =1 


1=1 


n—1 


k =0 

72 — 1 


fc +1 


+ 2 ^ a 2(n-fc-i) ( ^ 4- 2 )a 2 ja 2 ( fc _ z+2 ) 


z=i 


n—k 


n—k 


(3.7) 


-^(2fc + l)(2fc + 2)a 2fc+ 2(^ a 2l a 2(n-k-l)J + 7 ^ a 2 k{^22 a 2l a 2(n-k-l) y 
k= 0 1=0 k= 0 Z=0 

for all integers n > 1. The relations (3.6) and (3.7) show that a 2n +2 is uniquely determined by 
ao,, Qq , 7 ,..., 7 ( 2ti ) . This completes the proof. □ 
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